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ABSTRACT 

We evaluate the contribution of cosmic microwave background (CMB) polarization spectra to cosmological 
parameter constraints. We produce cosmological parameters using high-quality CMB polarization data from 
the ground-based QUaD experiment and demonstrate for the majority of parameters that there is significant 
improvement on the constraints obtained from satellite CMB polarization data. We split a multi-experiment 
CMB data set into temperature and polarization subsets and show that the best-fit confidence regions for the 
ACDM six-parameter cosmological model are consistent with each other, and that polarization data reduces the 
confidence regions on all parameters. We provide the best limits on parameters from QUaD EE/BB polarization 
data and we find best-fit parameters from the multi-experiment CMB data set using the optimal pivot scale of 
k p = 0.013 Mpc" 1 to be {h 2 n c , h 2 {l b , H , A s , n s , r}={0.1 13, 0.0224, 70.6, 2.29x 10" 9 , 0.960, 0.086}. 

Subject headings: CMB, anisotropy, polarization, cosmology, cosmological parameters 



1. INTRODUCTION 

The cosmic microwave background radiation contains fluc- 
tuations in temperature and polarization which have specific 
spectral features that record the evolution and constituent 
properties of the universe. The radiation is predicted to be po- 
larized at the 10% level due to Thomson scattering in the pres- 
ence of velocity inhomogeneities in the photon-baryon fluid at 
last scattering. 

The standard cosmological model predicts acoustic peaks in 
CMB intensity and polarization spectra. Polarized CMB ra- 
diation can be decomposed into even-parity E-modes which 
are generated by scalar and tensor perturbations, and odd- 
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parity B-modes which are generated by gravitational waves 
and gravitational lensing effects. In this paper we exam- 
ine the cosmological implications polarization spectra in- 
cluding the spectra produced by the 2009 improved anal- 
ysis of QUaD second and t hird season observations pre- 
sented in iBrown et alJ ((2009). QUaD i mproves on the de - 
tections of E modes made by the DASI dKovac et al.lf 2002). 
CAPM AP dHedman et alJ 120021) . Boomerang dJonesetal] 
120061). Wilkinson M icrowave A nisotropy Probe (WM AP) 
dDunklev et alJ 120091) . and CBI dReadhead et all 120041) ex- 
periments, adding accurate polarization data at small angular 
scales. 

We presented design and optics repor ts desc ribing 
the QUaD 1 experim ent in iHinderks et al.l (l2009h and 
IP' Sullivan e"tal]d2008l). Th e first season of QUaD results ap- 
peared in lAde et alJ d2008l) . We first reported measurements 
of temperature and p olarization from th e second and third sea- 
son QUaD data in iPrvke et alJ d2009l) and cosm ological pa- 
ramete r analysis of the data was carried out in ICastro et alJ 
(2009). The second and third season data sets from the QUaD 
experiment have been reanalyzed using a ground template 
removal technique rather than field-differencing to remove 
ground contamination. A description of the method and maps 
and spectra from two parallel a nd independent pipelines can 
be found in IBrown et alJ d2009l) . The re-analysis has effec- 
tively doubled the size of the QUaD field, increasing the pre- 
cision of the CMB spectra measurements by ~ 30% and con- 
straining the amplitud e of the lensing B- modes to <0.57 /j,K 2 
at 95% confidence. In lBrown et alJ d2009l) we also presented 
a cosmological parameter analysis which gauged the effect of 
using the improved spectra from QUaD in combination with 
WMAP and ACBAR data, demonstrating that QUaD data 
adds significant power for constraining cosmological models 
which include tensors and running of the spectral index. 

1 QUaD stands for "QUEST and DASI". In turn, QUEST is "Q & U Ex- 
tragalactic Survey Telescope" and DASI stands for "Degree Angular Scale 
Interferometer". 
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TABLE 1 

Baseline and Derived Parameters and Flat Priors on Baseline 
Parameter Set 



L/CALl 11J L1VJ11 


Pamm ftpr 

L 111 til 1 Id. CI 


Prior 


Baryon density 


n b h 2 


0.001-0.999 


Cold dark matter density 


Q c h 2 


0.001-0.999 


Acoustic peak scale 


e 


0.3-12 


Scalar fluctuation amplitude 


ln(10 10 A s ) 


2.7-4 


Scalar fluctuation index 


n s 


0.01-2 


Optical depth 




0.01-0.8 


Age 


Age (GYr) 


10-20 


Dark energy density 


n A 




Matter density 






Reionization depth 


Zre 




Hubble constant 


H 


40 - 100 


Linear mass perturbation 


08 








FIG. 1 . — QUaD polarization spectra (black) plotted together with WMAP, 
Boomerang, CBI and BICEP polarization spectra. The black lines are the 
QUaD+CMBmany theoretical best-fit spectra. 

In this paper we carry out cosmological parameter estima- 
tion in which we explore separately the strengths of CMB 
temperature data and CMB po larization data in cluding the 2 
year data release by BICEP (Chian get al.ll2.01QI) . highlighting 
the contributions of each to a multi-experiment analysis. 

2. DATA AND METHOD 

The QUaD data set consists of 143 selected days of data 
measured during the Austral winters of 2006 and 2007, at 100 
GHz and 150 GHz, in a region of approximately 100 deg 2 , 



and has produced polarization spectra of unprecedented qual- 
ity dPrvke et alj|2009t iBrown et al.ll2009l) . The QUaD experi- 
ment is the first experiment with the sensitivity to detect mul- 
tiple acoustic oscillations in the E-mode spectrum and TE- 
spectrum up to £=2000 in additio n to providing the lowest 
upper-limits on B-mode detection (Bro wn et al.ll2009l) . The 
spectra we used for cosmological data analysis, TT, TE, EE 
and BB, are an optimally weighted combination of three sets 
of spectra: 100 GHz, 150 GHz and 100GHz-150GHz cross, 
each with 23 band power values in the range 200 < I < 2000. 

Our CosmoMC-based parameter analysis is br oadly simi- 
lar to t he Monte-Carlo Markov chain analysis of ICastro et al.l 
(2009), hereafter Parameter Paper 1. A single, constant co- 
variance matrix was estimated from simulations within the 
QUaD pipeline. The covariance matrix is populated only in 
the diagonal, 1st off-diagonal and 2nd off-diagonal terms of 
each of the sub-blocks of the spectra. TT-TE and TE-EE have 
no 2nd off-diagonal terms and TT-EE and EE-BB have co- 
variance terms up to the 12th and 4th band powers respec- 
tively. Beyond these regions the covariance matrix would be 
dominated by noise from the numerical simulations. We also 
generated offset-lognormal x-factors from noise-only simu- 
lations. This enables us to model the offset-logn ormal likeli- 
hood f or our parameter estimation as suggested by lBond et al.1 
(2000). 

We make use of templ ates of the Sunyaev-Zer dovich 
(SZ) amplitude derived from Kom atsu & Selj ak (2002) which 
model a frequency-dependent contribution to the temperature 
power spectrum from the thermal SZ effect 2 . We apply a cut 
of spectral power above £=2000 to avoid the effects of resid- 
ual point-source contamination, and analytically marginalize 
over the SZ amplitude 3 . The justification of a fit to SZ pa- 
rameters on scales below £=2000 has not been confirmed 
by the high multipole temperature data obtained by QUaD 
dFriedman et alJl2009l) . 

We estimate likelihoods from our CMB distribution us- 
ing the publicly available CosmoMC dLewis & Bridldl2002l) 
Monte Carlo Markov chain algorithm. The theoretical CMB 
model at every stage is o btained from the p ublicly available 
CAMB Boltzmann code dLewis et aU HOQO). The choice of 
baseline parameter set, and the shape of the priors on param- 
eters fitted concurrently will impact on the one-dimensional 
(ID) marginalized parameters calculated using Monte Carlo 

2 SZ templates are available at http://lambda.gsfc.nasa.gov 

3 Our analysis in Parameter Paper 1 did not include marginalization over 
an SZ amplitude or offset-lognormal factors. 
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Markov chain. For the standard ACDM six-parameter cos- 
mological model, we limit this effect by using the base- 
line parameter combination comprising the baryon density, 
cold dark matter density, acoustic peak scale which records 
the expansion history of the universe, scalar spectral ampli- 
tude, scalar spectral index and optical depth: {/i 2 Sl/,,/i 2 f2 c ,#, 
ln(10 10 A. s ), Mj,t}, the use of ln(10 10 A. s ) also considerably in- 
creases the rate of convergence; Hq = 100/ikm s" 1 Mpc" 1 is 
the Hub ble constant. This is th e default parameter set of Cos- 
moMC (iLewis & Bridlell2002l) . 

When presenting two-dimensional (2D) plots of likelihood 
space evaluated from a CMB spectrum that does not have in- 
formation in the angular range below the first acoustic peak, 
we use the combined A. s exp (-2 r) parameter, a measure of the 
overall amplitude of the fluctuations in the spectra optimally 
sampled in the analysis using the baseline parameter set. The 
pivot scale used to evaluate the amplitude and spectral index 
is k p = 0.05 Mpc" 1 in the case where we use QUaD data exclu- 
sively. The pivot scale used for multi-experiment CMB spec- 
tra analysis is k p = 0.013 Mpc " 1 in common wit h our cosmo- 
logical parameter results of iBrown et al.l (120091) . We impose 
flat priors on the baseline parameter set (Table [TJ, make the 
assumption that the universe is flat and include in all analyses 
the effects of weak gravitational lensing. 

In addition we make use of the CosmoMC utility to analyt- 
ically marginalize over "nuisance" parameters, which is im- 
pleme nted for beam and calibration uncertainties (Brid le et al.1 
2002). The effective beam size for the combined QUaD spec- 
tra is 4. 1 arcmin. The beam error is a function dominated by 
the sidelobes and varies with scale. The calibration uncer- 
tainty is 6.8% in power. 

We provide parameter constraints on these baseline parame- 
ters and additionally include plots and constraints in the more 
traditional format of {h 2 flb,h 2 fl c ,Ho,A s ,n s ,T}. 

We present 68% and 95% confidence regions for cosmo- 
logical parameters and 2D and ID marginalized plots using 
our independent statistics code on 10 or more Monte-Carlo 
Markov chains with at least 100,000 converged steps. Tests of 
con vergence were carried o ut using the Gelman-Rubin statis- 
tic dGelman & Rubinll 19921) . We also made use of the Getdist 
statistics package, which is bundled with CosmoMC, as a con- 
sistency check of our results. 

We carried out a parameter fit from the products of both 
QUaD data analysis pipelines to confirm consistency. 

2.1. Other Data 

We comb ine QUaD with th e following CMB data sets: 
ACBAR dReichardt et al.ll2009l) . including the 2% beam error 
on the 5 arcmin beam, a nd calibration error 4.6% in power; 
CBI (ISievers et al.l l2009|) with a stated ca libration error of 
2.6% in power; BICEP (IChiang et al.ll20Tol) , 5.6% calibration 
error and beam err or varying with scale; the WMAP 5 year re- 
lease dNolta et al.ll2009t) using the publicly available WMAP 
5 year likelihood software 4 . We describe the set of CMB data 
above, excluding QUaD, as CMBmany. In Figure Q] QUaD 
polarization spectra are plotted together with WMAP, CBI, 
Boomerang 5 and BICEP polarization data. 

3. RESULTS 
3.1. Parameter Constraints from QUaD data 

4 WMAP likelihood software is available at http://lambda.gsfc.nasa.gov 

5 QUaD shares its entire survey region with Boomerang, therefore we do 
not include Boomerang spectra in CMBmany. 



We present parameter constraints obtained from QUaD data 
in Table [2] The table presents cosmological parameters cal- 
culated from the QUaD TE/EE/BB spectra subset, QUaD 
EE/BB spectra , QUaD TT spectra and the full QUaD data 
set. The parameters obtained from the subsets are consistent. 
Cosmological parameter limits from the CMBmany data set 
are also presented in Table [2] 

In Figure |4] we display the parameter constraints and con- 
tours obtained from QUaD polarization data (TE/EE/BB) 
alone. The resultant confidence regions are considerably 
smaller than the corresponding confidence regions from 
WMAP 5 year TE data, over-plotted in Figure @] for all pa- 
rameters except A s and r. The parameter constraints arrived 
at using QUaD TE/EE/BB are within 95% agreement with the 
parameter constraints from the full WMAP data set and those 
from WMAP TE spectrum. We have carried out, in addition, 
a combined analysis using QUaD and BICEP polarization 
data, these latest-generation ground-based polarization exper- 
iments detecting polarization at both low and high-£ The 
preferred parameter regions for QUaD and BICEP TE/EE/BB 
data are also presented in Figure [4] We note that the centre of 
the preferred range for the scalar spectral index from QUaD 
TE/EE/BB data is ~ 0.77. We also note that this preferred 
range in n s concurs with the QUaD polarization analysis of 
Parameter Paper 1, which found n s = 0.766 ±0.152. 




FIG. 2. — Likelihood contours obtained from QUaD pure polarization data 
(EE/ BB) at 95% and 68% levels for the scalar spectral index and the baryon 
density, displaying the tail of the distribution to high values of QUaD 
EE/BB only constraints red/orange dotted; QUaD EE/BB constraints with a 
Gaussian prior on the Hubble constant blue/green dashed; constraints from 
WMAP 5 year polarization data black solid lines. 

The best-fit spectral index value obtained from the QUaD 
EE/BB spectra subset is 0.53. We show in Figure [2] that the 
likelihood surface is skewed over the range of scalar spectral 
index values. We over-plot the likelihood for WMAP TE to 
demonstrate the best-fit regions for the spectral index over- 
lap. To supplement our QUaD EE/BB only parameter fits we 
also carried out an analysis where we implement a Gaussian 
prior on the Hubble constant with the latest indep endent lim- 
its ob tained by the Hubble Space telescope (HST. lRiess et alj 
120091) 74.2 ± 3.6kms _I Mpc" f . The Hubble constant primarily 
impacts the width and separation of the spectral peaks, while 
the spectral index is evaluated from the tilt of the data. Apply- 
ing the information we have on the Hubble constant improves 
our parameter confidence ranges, while not preferentially im- 
pacting the spectral index signal in the data. These contours 
are also plotted in Figure [2] confirming the robustness of 
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TABLE 2 

COSMOLOGICAL PARAMETER CONSTRAINTS USING QUAD DATA 



Parameter QUaD TE/EE/BB QUaD EE/BB QUaD TT QUaD CMBmany 11 

Thh 2 0.0233 +00030 0.0327 +0-00" 0.0218 +0.0040 0.0243 +°, gog 0.0225 +om( > 

n c h 2 0.124+S 0.162+™? 0.117+™! 0.119+°$$ 0.114 



+fW +8W S Mi° JoW 5 u.uzzj lqqq^ 

-0.030 0.162 _nmi 0.117 _nW 0. 1 19 -oioa 

1.040 ± 0.006 1-032 1.045 «[j 1.041 ^ lm1 h{Sh 

r b <0.53 (95 %cl) <0.46 (95 %cl) < 0.54 (95 % cl) < 0.54 (95 % cl) 0.087 +°-°}7 

ln(10 10 A t ) c 3.53 ;0 3| 3.74 +0f 3.48 +°-« 3.58+0-28 3.09 



i +0.117 n aru 



0.768 +«;!=o 0.530 0.920 0.804 +"S 0.964 1 013 

Age(GYr) 13.7 13.1 +«f 13.6 1°« 13.5 13.7^] 

f7 A 0.64±0.18 0.45+o| 0-68^ 0.68 +° 0.72 +°-* 

n,„ 0.36+°,!? 0.55+o : zl 0.32 +M 0.32+o : |3 0.28 +™ 



1+^.0 184+ 5 ^ 24 3+^°'^ 23 9+ 8 - 1 105+'- 4 

-10.0 10 - -9.8 zqo -11.2 ZJy -9.0 lu ">-1.4 

+12.5 « -1+13.1 70 s +15.8 71 1+10.9 7 n a +2.4 



Zre 23 

ffo 68.6™ 63.2™-| 72.5+f| 71.1 ™| 70.4 

CT8 d o.98±o.i7 i-07 +8:is o- 98 ^:!! o- 98 ^ °- 82 ^:o3 



a CMBmany describes the CMB data set including ACBAR (up to / =2000), CBIPol, BICEP and WMAP5 
(seeim. 

b The parameterization of reionization by CAMB has recently been changed. First ionization of helium and 
hydrogen are now modelled as simultaneous, affecting the returned values of r and the time at which z re is 
defined has changed at the 6% level. 

c The pivot point for A, and n s used with all the data sets presented in this table is k p = 0.05 Mpc~' . 

d The root mean squared linear mass perturbation is defined on a scale of 8/r'Mpc, where Hq = 100/ikm s 

Mpc~>. 



the EE/BB analysis. The parameter constraints from QUaD 
EE/BB spectra improve on the corresponding constraints we 
presented in Parameter Paper 1 . 

The constraint obtained on the acoustic peak scale, 9, is as 
tight from QUaD EE/BB only as it is from the QUaD temper- 
ature spectrum and the constraint on h 2 Q, c is tighter. 

The parameter ranges arrived at by using the complete 
QUaD data set agree with the values from the CMBmany data 
set within the 95% contour limits, as demonstrated in Fig- 
ure [5] The preferred parameter range for the bar yon density 
we pre sented in Parameter Paper 1 from the full iPryke et all 
(2009) QUaD data set was ~ 50% higher than that we now 
provide. The Hubble constant was also ~ 30% higher than 
the value we return from the analysis of the new non-field- 
dif ferenced data . The d ata analysis pipeline whic h produced 
the iBrown et alj d2009l) spectra differs from the Prv ke et all 
(2009) analysis in the implementation of a new noise-removal 
strategy and in modelling a new beam shape. This results in a 
material difference in the spectra in addition to a reduction of 
cosmic variance. The scalar spectral amplitude of Parameter 
Paper 1 (which we quoted as part of a combined A s exp (— 2r) 
parameter) was significantly shifted relative to our current 
preferred range, the corollary being in extended parameter 
analysis isocurvature ratios were more constrained. However, 
as our parameter analysis has evolved the confidence limits 
presented in Parameter Paper 1 cannot be directly compared 
to the limits in this paper. Repeating the para meter analy- 
sis of the second and third season QUaD data of IPryke et alJ 
(2009), using our current CosmoMC -based pipeline, we find 
68% confidence regions for A s exp(-2r) and the Hubble con- 
stant that are 30% greater than the c orresponding values ob- 
tained using the improved data set of lBrown et alJ (12009). 

3.2. Parameter Constraints Combining Data Sets. 

We continue our parameter analysis by investigating the ef- 
fect of adding QUaD data to the CMBmany temperature and 
polarization ensemble data set. 

In order to investigate the influence of polarization in 
our cosmological analysis, we carried out cosmological pa- 
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FIG. 3. — 2D marginalized contours h 2 £l/, against 9 from separate tem- 
perature and polarization combined QUaD and CMBmany data sets. QUaD 
polarization + CMBmany polarization constraints red/orange dotted; QUaD 
temperature + CMBmany temperature constraints blue/green dashed; con- 
straints from QUaD + CMBmany: black solid lines. 

rameter fits using the TT and TE/EE/BB subsets of the 
QUaD+CMBmany data set. The resulting best-fit regions, 
presented in Table [3] and Figure [3] are consistent, although 
the polarization spectra prefer lower values for the Hubble 
constant and acoustic peak scale as well as a lower spectral 
index. The multi-experiment TE/EE/BB data set contributes 
information on all six cosmological parameters. 



4. CONCLUSIONS 

We present best-fit regions for the standard ACDM six- 
parameter model using QUaD CMB spectra alone, and QUaD 
data in combination with CMBmany (BICEP, ACBAR, CBI 
and WMAP). 

QUaD data provide an unprecedented amount of indepen- 
dent information over four spectra in a multipole range of 
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TABLE 3 

Cosmological Parameter Constraints using combined CMB Temperature and 

Polarization 

Parameter QUaD + CMBmany TE/EE/BB a QUaD + CMBmany TT b QUaD + CMBmany 

0.0224 ±0.0005 
1 n ^ 005 

-0.QQ5 

1-04151 
0.086 « 7 

T 1 T +0.04 

-0.03 

0.960 
13.6*| 
73 +o iB 
0.27 
10.4 f{ 
70.6 +2:3 





0.0204 ±0.0020 


0.0231 ±0.0008 


Q c /i 2 


0.124 ±0.014 


n , no +0.007 

i n42 +0.003 
-Q.QQ2 
r2 7 +0.075 
u - li; -0.079 


e 


1.036 ±0.005 


T c 


n-i-i +0.017 
U.u/J _ on 


111(10'%) d 


o 97 +0.08 
J - zz -0.08 


3 17 ml5 

-0.15 


« s d 


0.842 ±0.055 


903 +0.026 


Age (GYr) 


14 1 -* - 3 

-0.3 


1 3 6 +0 - 2 


n A 


0.63 ±0.10 


75 + 04 

-QQ4 


n m 


0.37±0.10 


0.25 « 
12.6 


Zre 


10-1 


Ho 


63.8;|1 


73 7 +3.8 


CT8 


79 +fl 06 

u -' y -0.07 


0.83 |t 



0.82 



0,(13 
1.4 



-0.03 



a CMBmany TE/EE/BB is the combination of WMAP, CBI and BICEP polarization spectra. 

b CMBmany TT is the combination of WMAP, CBI, ACBAR and BICEP temperature spectra. 

c The parameterization of reionization by CAMB has recently been changed. First ionization of helium and 

hydrogen are now modelled as simultaneous, affecting the returned values of r and the time at which z r e is 

defined has changed at the 6% level. 

d The pivot point for A s and n s is k p = 0.013 Mpc~' . 



200 < I < 2000. This improves the ensemble of CMB data 
sets specifically providing independent cohesive information 
over an angular range broad enough to span the sparsely pop- 
ulated "hinge" multipole regions of the CMBmany data set. 

Data from the QUaD experiment, as a single, discerning 
data set can provide excellent independent measures of the 
baryon density in the universe, the cold dark matter density, 
and the acoustic scale, which measures the universe's evolu- 
tion history. 

QUaD TE/EE/BB spectra provide considerably tighter con- 
straints on four of the six standard ACDM parameters than 
can be obtained from 5 year WMAP polarization data. We 
present the best confidence ranges obtained from QUaD 
EE/BB spectra alone. The QUaD temperature spectrum alone 
also provides good parameter constraints in excellent agree- 
ment with the WMAP 5 year best-fit parameters. 
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FIG. 4. — Parameter constraints and contours obtained from QUaD polarization data (TE/ EE/ BB) and QUaD polarization combined with BICEP polarization 
using a pivot scale of k p = 0.05: (a)2D contours at 95% and 68% levels for the traditional physical parameter set making use of the combined A s exp— 2t 
distribution; (b) constraints on the six traditional physical parameters; (c)lD marginalized distributions of Monte Carlo Markov chain baseline parameters; QUaD 
polarization only constraints red/orange dotted; QUaD + BICEP polarization constraints blue/green dashed; constraints from WMAP 5 year polarization data 
black solid lines. 
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FIG. 5. — Parameter constraints and contours obtained from QUaD temperature and polarization data using a pivot scale of k p = 0.05: (a)2D contours at 
95% and 68% levels for the traditional physical parameter set making use of the combined A s exp-2r distribution; (b) constraints on the six traditional physical 
parameters; (c)lD marginalized distributions of Monte Carlo Markov chain baseline parameters; QUaD temperature and polarization constraints red/or ange 
dotted; QUaD + BICEP temperature and polarization constraints blue/green dashed; constraints from CMBmany: (ACBAR/ WMAP/ CBI/ BICEP) (see |2.it 
black solid lines. The best-fit value of r from QUaD temperature and polarization data is 0.314. 



